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Abstract 

The stationary states of boundary driven zero-range processes in random media with quenched 
disorder are examined, and the motion of a tagged particle is analyzed. For symmetric transition 
rates, also known as the random barrier model, the stationary state is found to be trivial in absence 
of boundary drive. Out of equilibrium, two further cases are distinguished according to the tail of 
the disorder distribution. For strong disorder, the fugacity profiles are found to be governed by 
the paths of normalized a-stable subordinators. The expectations of integrated functions of the 
tagged particle position are calculated for three types of routes. 
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I. INTRODUCTION 

The zero-range process, first introduced by Spitzer [l!], is a simple model for a gas of 
mutually interacting particles on a finite or infinite lattice. The interaction is local in the 
most strict sense: the rate, at which particles leave a lattice site, depends on the number 
of particles present on that particular site only. This restriction makes the stationary state 
analytically tractable-it is a product measure-even in presence of driving fields that break 
the spatial symmetry of the transition rates. The invariant measures of the process on TL^ 
were studied by Andjel {2]. We shall consider finite lattices only. 

Due to its analytical tractability, the zero-range process has been used to test ideas con- 
cerning the stationary states and dynamics of real interacting systems. The hydrodynamic 
limits were given a rigorous treatment in the book of Kipnis and Landim [sl, which provided a 
basis for recent studies concerning the thermodynamic functionals of nonequilibrium steady 
states ^. The applications of zero-range processes in physics range from shaken granular 
gases to condensation phenomena (see [sl and references therein). 

The stationary states of zero-range processes with open boundaries have already been 
discussed in detail by Levine, Mukamel, and Schiitz jo]. In this article, we shall consider 
open systems but the media itself will be a random in the sense that the time evolution of 
the particle system takes place in a fixed environment, but the exact structure of the media 
is not known-only its statistics for a generic sample. This type of randomness is known 
as quenched disorder. Bulk driven zero-range processes in random media have been used 
io study jamming transitions in asymmetric exclusion processes with particle- wise disorder 
1, Q]. We shall show that there is plenty to discover even in the boundary driven but 







otherwise symmetric system. In particular, the stationary state shows nontrivial correlations 
when the system is driven out of equilibrium. By restricting the dynamics to the most 
simple and natural case of symmetric media, where the jump rate from lattice site i to i + 1 
equals the rate in the opposite direction (imagine the particles surmounting barriers with 
independent and identically distributed random heights), we shall find that, due to universal 
laws for normalized sums of independent and identically distributed random variables, closed 
expressions for a number of quantities can be found. 

In addition to characterizing the structure of the stationary states in random media, we 
shall also focus on the motion of a single, tagged particle in a gas of zero-range interacting 
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particles. This problem is on one hand connected to the theory of diffusions and random 
walks in random environments reviewed in j^, [lo|, 11 1 (Isichenko 1^ gives emphasis to the 
motion of a passive tracer in random flows). Most of the studies concerning random walks in 
random environments allow the asymmetry of the transition probabilities, which makes the 
problem much harder for zero-range processes. Even the model with particles diffusing in a 
landscape of valleys with random depths (the particle jumps out of a site to either direction 
with equal probability [9]) presents some extra difficulty as compared to our model with 
barriers of random height, because the existence of the stationary state is, in that case, not 
easily controlled. Rigorous results for a single particle in symmetric random environment 
were derived by Kawazu and Kesten |l^ • They proved that the suitably speeded up random 
walk converges on the level of path measures to an a-stable diffusion as the process is 
observed on larger and larger length scales. Such convergence results are sometimes called 
"invariance principles" |l^ . 

On the other hand, there exists a vast literature on limits of tagged particle rnotion in 
exclusion and zero-range processes, and some other interacting particle systems [14|, mostly 
in homogeneous media. There are a few important contributions for zero-range processes: 
For symmetric and translation invariant transition probabilities, Saada [15] proved a cen- 
tral limit theorem for the position of the tagged particle for a particular interaction, and 
Siri [l^ showed an invariance principle for processes on tori, with the interaction possibly 



depending on the location. Sethuraman 171] considered the variance of the particle position 
in asymmetric cases. 

We shall show how the expectations of certain functionals of the paths of the tagged 
particle, more precisely functions of the tagged particle position integrated up to stopping 
times of an auxiliary process, can be calculated first in the finite system and, in the end, as 
the macroscopic limit is approached. The convergence of the motion in random symmetric 
media to a diffusion on a formal level is not proved. 

The structure of the article is as follows. In the next section, we define the model and 
solve its stationary distribution in a fixed environment. For later use in this article, we 
also derive the generator of the adjoint dynamics. Section IllBI discusses the stationary 
states in random media. The cases of weakly (III B 1|) and strongly (III B 21) inhomogeneous 
environments are treated separately. Some properties of the fugacity profiles in the latter 
case are discussed in section 111 B 31 The section Hn] is devoted to the motion of a tagged 
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particle, again first in a given environment (section IlII Ap . Its movement in random media 
is considered in IIIIBI The results are commented and related to other studies in the last 
section. 

Throughout the article, the symbols a V 6 and a A 6 are used, respectively, to denote the 
maximum and the minimum of two numbers a and 6, and [xj stands for the greatest integer 
smaller than or equal to x. 



II. STATIONARY DISTRIBUTION 



A. Preliminaries 



The nearest neighbour zero-range process {X{t))t>Q = (Xi(t), . . . , X„(t))(>o on [n] = 
{1,2,..., n} with left (resp. right) boundary reservoir at (resp. + 1) is a Markov process 
with values in Z" for each t > 0, and generated by 

n 

+Y.9iv^){v^[fW''^') - m]+umv'''-') - fiv)]}, (1) 

i=l 

where rj'^'^ is the state obtained from r] by moving a particle from i to j with the convention 
that ?7°'^ = (?7i + 1,T]2, . . . ,T]n) etc. for transitions involving boundary reservoirs. The single 
particle rates {ui, Vi) for the jumps to the left and right, respectively, are assumed positive and 
bounded, and the interaction g : Z+ — > IR+ is such that, for some e > 0, = (?(0) < e < g{k) 
for A; > 0, and has bounded increments; sup^ \g{k + 1) — g{}i)\ < K. 

The next proposition involves the mean current from j to j — 1, which is defined as the 
expected number of particles that move from j to j — 1 during one unit of time minus the 
expected number of particles that move in the opposite direction. We shall also use the 
shorthand notation g\{m) = YYjLi qU) for m > and (71(0) = 1. 

Proposition 1. The stationary distribution of the zero-range process generated by ([1]) is 
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the product measure fi{r]) = YYj=i ^<l>jiVj) with 

0. = ^ri^+-i:ri^, (3) 



^ ' Z^i=o llfc=o «„_fc 

given that the fugacities satisfy (pj < hminfm^oo for j G [n]. Here c equals the 
mean current from j to j — 1 for any pair of sites with j G [n + 1] . 

Proof. We show that -E'^[(L„/)(X(t))] = = for bounded /. By a change of sum- 

mation index, this holds if /i satisfies the balance equations 

n 

for T] G Z". Substitution of the product form with ([2]) and some simple algebra leads to 
traffic equations for the fugacities 0^: 

- Vj(f)j = Uj(j)j - Vj^i(f)j_i (6) 

Ul4>l -Vq = Un+l - Vn4>n, (7) 

the solution to which is given by and (jl]) with c equal to the quantities on both sides of 
the equations. The expected number of particles that leave site j in a time interval of unit 
length is J^rii'^j + '^j)9i''lj)f^{v) = (^j + '^j)4>j- So c really is the mean current. □ 

Given that (pj < liminfm->oo for j G [n], the expected number of particles at site j 
is finite and given by the expectation of a i/^^ -distributed random variable, 

p, := R{<f>,) := <I>,Z' /Z (8) 

Notice that the function R is strictly increasing because (f)R'{(f)) is the variance of the distri- 
bution z/0, which is positive by the assumption that the increments of the interaction g are 
bounded. 

Formula (jl]) shows that, in general, there is a current of particles due to the imbalance of 
the fields (uj) and (wj), which means that the process is asymmetric with respect to reversal 



of time. Since the time reversed process is generated by the adjoint of L„ in L'^di), that is 
ljL[h{Lnf)] — //[(-L*/i)/], change-of- variable calculations lead to the following: 

Proposition 2. The time reversal of {X{t)) with respect to is a zero-range process with 
the generator 

Kfiv) = uMfiv''') - fiv)] + VnUfir^''^ - fiv)] 

i.e. it has single particle rates v* — and u* — with the convention 

that 00 = = 1- 

This shows that the process is reversible if and only if the detailed balance conditions 
Vi(f)i = Ui+i(f)i+i hold, i.e. c = 0. 

In this article, we shall restrict our attention to symmetric media. This means Uj+i — vj 
for j' e [n — 1] and possible exceptions on the boundaries, where we take Vq = v as the input 
rate to the system from the left reservoir and Un+i = u from the right. The output rates at 
the boundaries are chosen to equal unity, ui = Vn = 1- Define 

«. = i + E^r', (10) 

1=1 

with si = 1. Then proposition 1 yields the following: 

Corollary 1. For symmetric media with boundary fields u, v, the stationary distribution is 
of the product form as in Proposition 1 with 

c|>^ = U^+v(l-^), (11) 



U — V 



given that 0i V 0^ < liminf g{k). 



(12) 



Naturally, the process is time reversible if and only if the boundary rates coincide, and 
then the fugacity takes a constant value equal to the input rates. For general u and v, (pi is 
a monotone function of the spatial position, and therefore a condition only at its boundary 
values need be imposed to guarantee the existence of the stationary distribution. In the 
following sections, we shall adopt a stricter condition, u\/v< liminf g{k), than the one used 
in the corollary. 
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B. Symmetric random media 



We shall now consider zero-range processes in symmetric random environments with 
quenched disorder. This means that the environment does not evolve in time. For this 
purpose, let the single particle rates for jumps between i and z + 1 in the bulk be 

independent and identically distributed random variables on some probability space with 
measure P. Expectations with respect to this measure will be denoted by E. The symbols P 
and E are reserved for the Markov process in a fixed environment. Again, the assumption of 
positive and bounded rates is made. Each random environment is equipped with the same 
boundary rates as in the setting of corollary 1: The input rate Vq from the left reservoir to 
site 1 is denoted more conveniently by v, the rate Un+i from the right reservoir to site n by 
u, and the output rates from sites 1 and n to the left and right reservoirs, that is ui and Vn 
respectively, are set to unity. The sufficient condition 

u \/ V < limini g{k) (13) 

for the existence of the stationary state is always assumed, and unless stated otherwise, the 
equilibrium case u = v is excluded from considerations. 

1. Weak disorder: E,v^^ < oo 

The existence of expectations of the inverse bulk rates determines the large scale structure 
of the system under study. In fact, the situation is rather simple for Ew^^^ < oo: Defining a 
rescaled continuum analogue of Eqn. flTUl) by 

L(n+l)xJ-l 
i=l 

and choosing a = 1 yields that Sn.ii^) xKv~^ for P-almost all environments as n — >■ 
oo. The convergence is even uniform because the functions Sn,a{x) are increasing, and the 
limiting function is continuous and increasing. Thus the continuum version of the fugacity 
satisfies 

uniformly almost surely. Moreover, the current is inversely proportional to the system size, 



The continuity of the map R{(f)) = 4>Z'{(f))/Z{(f)) on [u Av,u\/ v] imphes the uniform conver- 
gence of the particle density, 

p„(x) := {R o 0„)(x) — > R{ux + v{l — x)) (17) 

P-almost surely. As shown in [4j], the fluctuations of the macroscopic density profile are 
described by the free energy of an equilibrium system with spatially varying mean waiting 
times E(fj + Ui)~^ used to create the same fugacity profile. This should be contrasted with 
more complicated non-equilibrium systems, such as exclusion processes, with non-local large 
deviation functions llSl. 



2. Strong disorder: Ew^ 



oo 



We now turn to the nontrivial case of Ef ~^ = oo. In particular, we assume that (f 
belong to the domain of attraction of a strictly a-stable distribution |l9l . l20l| . with a G (0, 1), 
and 'strictly' meaning absence of a deterministic component as in Sato [19| . In other words, 
there is a slowly varying function L such that (nL(n))^^^°' Yl^=i '^1^ converges in distribution 
to a strictly a-stable random variable S* as n ^ oo, or equivalently (see [2^ for a proof), 
the tail of the distribution of the is heavy with exponent a. For simplicity, we take 
L{n) = 1. The sums of v^^ being almost surely positive, the properly rescaled current 
converges in distribution by theorem 5.2 of reference 21| , 

n^/^c — >{u- v)/S. (18) 

There is also a limit theorem for the fugacity as a stochastic process. The processes 
{Sn,a{x),x E [0,1]}, defined by equation (HM . as elements of the space -D[0, 1] of right- 
continuous functions on [0, 1] with left-hand limits, converge in the Skorohod topol- 
ogy (implied by the metric d((t)_, ip) _ = iiif{sup^ \(t>(x) — iIj{\{x))\ V sup^ |a; — A(a;)| 
A continuous and one-to-one) 
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22 



23 



creasing Levy process 
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24( 1) to a strictly a-stable subordinator (in- 



25|) {S{x),x G [0,1]}. The set of discontinuities of the map 



{Sn,a{x),x G [0, 1]} 1-^ {Sn,a{x) / Sn,a{^) , X G [0, 1]} from the spacc D[0, 1] to itself consists, 
in this topology, only of the path that is identically zero, and is thus of null measure for the 
process 5*. Then, by theorem 5.1 of reference [211], the fugacity processes (pnix) converge in 



distribution to 



(f){x) 



5(1) 



(19) 



where the normahzed strictly a-stable subordinator S{x)/S{l) on [0,1] is also a random 

known in mathematical statistics ClS db part 



26| . The next theorem gives the one-point 



distribution function. This object is rather wel 
of a family generalizing the Dirichlet process 
statistics of the fugacity. 

Theorem 1. The distribution function of the fugacity at x G (0, 1) is 

^, , , , ^, 1 sgn(u — v) / x\u — 9\" — (1 — x)\9 — vl"" 7ia\ 

FU(x) <e) = ^ ^ arctan ^ ^ ^- ^ tan — 20 

and its moments can be calculated from the formula 

n X \ k 

E<p{xr = ^" + 5^ ( ~ v)'Y,^'^c.,o{k,j), (21) 

k=l ^ ^ i=l 

where 

^ - i):^-Eri "^-°"^:!,nlr^-°»'°' - (-) 

i=l * V y 

The summation is over the sequences of nonnegative integers {mi)^^-^ such that Yli=i ~ J 
and Yl^=i ~ ^• 

Proof. Since the normalized subordinator 5'(x) / 5'(1) is a random distribution function on the 
set of real numbers, there is a random probability measure, to be denoted by ip, associated 
with it. Regazzini, Lijoi and Priinster '27| have derived expressions for distribution functions 
of integrals with respect to such measures using the Gurland inversion formula 

1 ^e^'^ 
lim 

for a distribution function F and the corresponding charateristic function (p. Since for u > v, 



F{0 + F(e-) = 1 - 1 hm ^{t) dt (23) 



/ N ^ S(X) 9 — V n 9 — V 



choosing / = /{■ < x} in formula (10) on page 575 of reference [271], which in our case reads 
^[P(^[/]<^)+P(^[/]<^^)] = i 

- ^ arctan { M) - Omv) - e\^^y ^^ 

and integrating yields the first statement of the theorem. Notice that there is null probability 
that (j){x) takes a fixed value for x e (0, 1), and thus ¥{(p{x) < 6) = P(0(x) < 6). 



To prove the second statement, by the binomial theorem we need to show that the k-th 
moment of the normahzed subordinator at x equals the generating function for the sequence 
Q'a,o(k,j), which is in fact a probability distribution, as will be explained in remarks after 
the proof. The use of double subscript shall also be justified. The statement is verified by 
the following calculation based on conditioning, Fubini theorem and Faa di Bruno formula. 
Here {J-'x)x>o is the natural filtration of the subordinator S{x), and $(A) = A"<l'(l) is the 
Laplace exponent of its (infinitely divisible) distribution at x = 1. 



E[S{x)''S 
= E 



E 



S{x)''E 



e-^«^^^MAi---dA, 





(-1)' 
(-1)' 



■ ■ ■ io ( gAi- ■c}a/ "^''^'''''0 e"^'"^^*^''^''^ dAi ■ ■ ■ dA. 



(9Ai ■ ■ ■ dXk 



{k-l)\ 



Q (A;-l)!Vdr'= 

^-ig-r-"* (1) ^ 



k 

mi - ■ -ruk 

k 



i=l 



dr 



1=1 



where the summation without indices on the two lines before the last one is over the set 



{("^i)f=i : "T-i > 0, Ylii=i i'lTT'i = ^} and = a{a - 1) ■ ■ ■ (a - z + 1). 



□ 



The double index notation in ga o comes from the theory of exchangeable random parti 



tions, and is taken from J. Pitman's lecture notes 



topic include 



29 



30 



31 



26|. The list of research articles on this 



32l |. We quickly review here the meaning of the sequence qafi{k,j) 



in the context of partitions. 

An exchangeable random partition is a partition of the set of positive integers N generated 

by sampling from a random discrete probability distribution. A partition is identified with 

the sequence of partitions of truncated sets [n] = {1, . . . ,n} in such a way that the probability 

of a partition of [n] in k parts with cardinalities {ni, . . . ,nk) is given by the symmetric 
! ' 

function [30] p{ni, . . . , n^) = jf.) ^ Y\i=i ^jl*- Here the summation is over all sequences 

of k distinct positive integers, and Pi is the i-th largest atom of the random sampling 
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distribution. One way of thinking here is through the species samphng models of ecology: 
A sample of n individuals from a large population, with ranked relative abundances given 
by the sequence (-Pi), is partitioned according to the species in that sample. One possibility 
is to take (Pj) from the statistics of jumps of an a-stable subordinator, which provides the 
connection to the system under study. A more general scheme known as the two-parameter 



family of Poisson-Dirichlet distributions, and denoted by PD(a, 9), is presented in j3l|. The 
number qa,e{k,j) is then the probability that in a sample of size k from a population with 
frequencies PD{a,9), there are exactly j species. 



Poisson-Dirichlet distributions appear naturally in the theory of spin-glasses 33|]. In 



particular, the ranked Gibbs weights have this distribution in the thermodynamic limit of 



Derrida's random energy model 3J]. PD(l/2,0) also appears as the statistics of ranked 
lengths of Brownian excursions, and PD{a, 0) for excursions of more general objects such as 
Bessel processes. Therefore one expects a wide range of applications for these structures in 
the physics of random systems. 

Example. As a specific application of theorem 1, we now determine the density profile near 
the boundaries for the attractive interaction g{k) = I{k > 0}. For u,v < 1, we get 
= 0(x)/(l — 0(x)) = By theorem 1, and the fact that the expected 

number of species in a sample of size k with PD(a;, 0) frequencies is (a -|- l)(a -|- 2) ■ ■ • (a -|- 



A;- 1)/(A;- 1)! |26|, 



k=l ^ ^ j=l 



k k 

JQafiikJ) 



u — V 



{l-u){l-v) \l-u^ 
On the other hand, the derivative at the opposite boundary, x ^ 0, is 

because qa,o{k, 1) = (1 — a)(2 — «) ■ ■ • (A; — 1 — a)/(A; — 1)! by formula (!22l) of theorem 1. 



3. Properties of the fugacity distribution 



Let us now discuss some consequences of theorem 1. First of all, the result for the 
moments shows that the expected value for the fugacity at x is K(j){x) = v{l — x) + ux. The 
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formula E0(x)2 - (E0(x)) = (1 - a){u — f — x) for the variance allows us to estimate 
the magnitude of deviations around the linear profile; 



P 



— t>(l — x) — ux\ > |n — v\^/x{l — x)^ < 1 — a. (26) 



This matches with the intuitive picture of the fugacity process: By the Levy-Ito decomposi- 
tion, the discontinuities are countable and dense on the spatial axis, but for a close to 1, the 
jumps do not differ in size very much, and the profile is usually nearly linear. For a tending 
to zero, a small number of jumps dominates the normalized sum in the fugacity, and the 
system is segregated into regions of nearly constant particle densities separated by bonds of 
very low conductivity. 

More information on the statistics of the fugacity can be extracted by differentiating the 
distribution function fl20|) . For convenience, we consider the case v < u only, and introduce 
the notation X = 6 — v,C, = u — d. Then the probability density reads 

f (0) = ^^^^^ (m - v)x{l - x)A"-^r~^ .271 

■''^^"^^ ' TT (xA" + (l-x)e)'cos2^ + (xA"-(l-x)e)'sin2ni" ^ ' 

It resembles the generalized arcsine distribution that appears as a solution to some occu- 
pation time problems. A common feature of the two densities is that they both diverge 
at the boundaries, which is perhaps unexpected for a close to 1 in the present case. The 
main difference is that the arcsine density is convex, whereas the density given above has 
two minima for large enough a. Taking into account the divergence at the boundaries and 
confinement within ux + v{l — x) ±{u — f)A/x(r^-^ as a — > 1, this is no longer a surprise. 

It is easily seen that the number of extrema of the density (1271) does not depend on 
the actual values of the boundary fields, as long as they are not equal. So without loss of 
generality, we choose v = and u = 1. The number of extrema within (0, 1) is then given 
by the number of solutions to 

(2^-l-a)02"(l-a;)2 + (2e-l + a)(l-0)2V + 2(2e-l)(l-0)°r(l-x)a;cos7ra = 0. (28) 

Let us concentrate for simplicity on the statistics in the middle of the system, i.e. x = 1/2. 
Writing 5 = 6^ — 1/2, we get 

^ + S^ (25-a)+Q-(5^ (25 + «) + 4(5 Q + (5^ ~ ^°^™ = °' ^^9) 

which has a solution 6 = for every a G (0,1). Moreover, the left-hand side is an odd 
function of 5, so if nontrivial solutions exist, they are found symmetrically from both sides 
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of the origin. Notice also that the derivative of the density is negative for S just above —1/2, 
and positive for 6 close to 1/2. There are at most three solutions. Let a* ~ 0.5946 be the 
solution to 1 — 2q;^ + costtq; = 0. Then the expansion of the left-hand side of equation (12^ 
to the linear order in 6 yields that for a < a* the density is convex, and a local maximum 
exists at 5 = in the opposite case. 

As is already clear, the particle densities measured simultaneously at distinct points are 
correlated because the structure of the underlying media is unknown. The correlations in 
the fugacity profile reflect this fact. In principle, it is possible to evaluate general multiple 
point correlation functions for the fugacities using the expansion 

n n k Q/ \ 

1=1 k=0 j=l ^ ' 

where the summation is over all subsets of {xj} of cardinality k. The calculations, using the 
same techniques as in the proof of the theorem above, are rather tedious for general n, and 
we just mention the result for the covariance. For < a; < ?/ < 1, 

E0(2;)0(?/) - E0(a;)E0(?/) = (1 - a){u - vfx{l - y), (31) 

which shows positive correlations. 



III. MOTION OF A TAGGED PARTICLE 



We shall first study the motion of a tagged particle in a boundary driven gas that has 
reached its steady state in fixed environment, but eventually we shall consider the case 
of random media. The convergence of the tagged particle processes to one- dimensional 
diffusions characterized by scale functions and speed measures (as explained for example in 



35, 



36| ) is not within our reach, but we content ourselves with the analysis of expectations 



of functionals of the particle position integrated up to the time of exit from a set. 



A. Travels in fixed environments 



Let the symmetric environment be fixed, and let J = {a + 1, a + 2, . . . , 6 — 1} C [n]. The 
particle system (Xi)^^^ is assumed to be in the stationary state at t = 0. Three types of 
paths of the tagged particle corresponding to three experimental situations are considered: 

13 



In the first case, a transition of tlie particle from a to a + 1 is observed at t = 0, and some 
time later, its first exit from J by a jump from 6 — 1 to 6 is registered. It is observed to 
return to a without reaching h in the second case. Finally, the starting position i G / is 
known in the third setting, and the time of exit from J by a jump through either end is 
measured. 

The position of the tagged particle will be denoted by Yt and the exit time from J is 
Tj = inf{t > : Yt = a or b}. Kook and Serfozo [sT] have constructed a process such that 
the travel times of the preceding paragraph are stopping times. This process keeps track 
of the paths that the particles are traveling, and labels them accordingly. Let P"-^^^ pa^b 
and P* be the Palm measures corresponding the three cases respectively (see 38| for Palm 
calculus). For example, 

P'^'iB) = \. / HY.^t e B)K^,{dt), (32) 

where Na^b counts the particles that, in the stationary state, start the trip from a to 6 and 
eventually complete it. 

It should be noted that the process of numbers of particles does not determine the motion 
of the tagged particle because it does not tell in which order the particles leave a site. In 
zero-range processes, it is usually implicitly assumed that the particle to leave the site is 
chosen in uniform random manner from the particles at the site just before the transition 
takes place. The results in this article hold for any rule that assigns a fraction fj, with 
X]j>i fj ~ 1' total "work" g{Xi) to the j-th particle of the total of Xj particles 

present, ranked by the order of their arrival to the site i. Notice that this extra freedom 
does not induce biased movement of a tagged particle. 

We now proceed to the calculation of expectations a class of functionals of the tagged 
particle paths under the Palm measures P"-^^^ pa^b^ ^^^^^ pi rjj^g functionals are of a 
specific type common in physics: They are functions of the particle position integrated up 
to the time of exit from a set under consideration. 

Define the continuum version of the function Sj in formula ffTOj) as s{x) = 1 + X]j=i^^ '^i^ ■ 
We shall also need the functions 

^ / X \s{x /\y) — s{a)]\s{h) — s{x\/ y)] , 

Gj{x, y) = slh) -s{a) ^ ^ ^ ^ ^ 
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and measures 

mn{dx)=J2^MWx), (34) 



i=i ^ 



where ((){ is the fugacity of equation (ITTll . Z is the partition function from equation ([2]), and 
6i is the Dirac measure at i. 

Theorem 2. For a function / : Z — > M, 

E'^-' / ■ fiY,)dt= / /(y)Gj(y,y)m„(dy) (35) 

Jo J-OQ 

E^-' r f{Yt) dt = r f{y) "if^'"^^^^] ^,y, rn^{dy) (36) 

Jo J-oo (s(6) - s(a))(s(6) - s(a + 1)) 

l-Tj POD 

E' f{Yt)dt= f{y)Gj{i,y)m^{dy). (37) 

Jo J-oo 

L 39I. 



Proof. Let us start with the first assertion. The proof is based on the Little laws 
for general travel times in queuing networks as expressed by Kook and Serfozo |37j]. Tj{i) 
denotes the time the particle spends at i before Tj, and pi{a h) is the probability that a 
particle at i is traveling from a to 6 without return to a. Then 

-^0 i=a+l 

= 15 fm'^-'Tjiz) = 15 f{{)pM b) J^^. (38) 

i=a+l i=a+l ^^^^ ' ' 



where the last equality holds due to theorem 3.1 in reference [37|. According to corollary 3.4 
in the same reference, the probability pi{a b) can be expressed as product ai{a, b)a*{a, h) 
of the probabilities that a particle at i is absorbed at h before hitting a, and that a particle 
obeying the adjoint dynamics (this is why we calculated L* in the preliminaries section) is 
absorbed at a before reaching h. Clearly, the first of these satisfies 

ai(a, h) = piai+i{a, h) + qiai_i{a, h) (39) 

with boundary conditions aa{a,h) = 0, ai,{a,b) = 1, and the shorthand pi = fj/(fi_i + fj), 
Qi = 1 —pi- By proposition 2, an equation similar to (!39ll holds for the probabilities a* with 
p* = V* /{y* + u*) = Vi(j)i+i/[{vi_i + Vi)(j)i] instead of Pi, and q* = 1 — p* instead of g^. The 
boundary conditions are replaced by a*{a,b) = 1, al{a,b) = 0. The solution to this set of 
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equations is 

a,(a,6) = ^^, (40) 

06 - (pa 
(piWb - (Pa) 

Noting that EXi = (j)iZ'{(j)i)/Z{(f)i), and ENa^ti^, 1] = "i^a^a^a+i some simple algebra 
finishes the proof of the first part. The second statement is shown in the same fashion. 
Observe that Pi{a ^ b) = ai{b, a)a*{a, b), and ENa^bi^, 1] = Va4>a(^a+i{b, a). 

The third part of the theorem can be shown to hold as a corollary to the other two. For 
a < i < j < b, 



E'Tj{j) = E'^'Tj{j) 



Vi-iai-i{i, a) + Viai+i{i, b) 



n>0 



Vi-i + Vi Vi-i +Vi J V Vi-i + Vi Vi-i + V, 



kj \Vi^iai_i{a,i) + Viai+i{b,i) J \Vi_iai_i{a,i) + Viai+i{b,i) 



X 

k=l 

Vi-iai-i{i, a) + Viai+i{z, b) Vi-iai-i[i, a) + Viai+i[i, b) 

In the first equality, the time spent in j while traveling from i to either boundary is split 
into the time in j for the trip from i to the right boundary without return, and into the 
time in j for the paths that do return to i. The outermost sum in the latter part is over the 
number of the returns. The inner sum then counts the expected number of right excursions. 
The left excursions do not contribute because i < j. The cases with j < i can be treated 
similarly. Substitution of the formulae for the expected occupation times and simple algebra 
now finishes the proof. □ 



B. Travels in random media 



As an application to the theorem of the preceding section, we now study the behavior of 
the tagged particle on macroscopic spatial intervals. To this end, we define J„ := { [an\ + 

1, . . . , [fenj — 1}, where < a < 6 < 1. The results of this section are presented informally 
and the proofs, or at least the sketches, are given within the text. 
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1. Weak disorder: Ev^ < oo 

Let us again deal with the case of finite expected inverse rates first. To get nontrivial 
hmits for the integrals in theorem 2 as n ^ oo, the motion of the tagged particle must 
be suitably speeded up. For this purpose, we shall consider the processes {Yn2t)t>o- Notice 
that the exit time T) of this speeded up process from a set J is identical in law with 
n~'^Tj. Let / : [0, 1] ^ M be continuous. Recall that by the assumption f|T3l) on the 
boundary rates, the measures m„ are bounded, with the weights of the Dirac measure 
given by a uniformly continuous function of the fugacity 0„(x), and that the functions 
Sn^i{x) = (1 + Yl\=i^^^^'^ '^^^M'^ converge uniformly almost surely. Then, by equation (I57|) 
of theorem 2, 

^^^^^ Jo " ^ U^"'7 ^ J ^^^Z'')^""^ La^^'J ' """^"^^^ 

^ ^ r/(y) ^^^"7^^'"^^"^ Mdy) (43) 

u-v Ja 0- a 

for almost all environments as n — oo, where m(d?/) = [Z'{uy + v{l — y))/Z{uy + v{l — 
y))] dy, because the integrand of the expression after the first equality, written in terms of 
the macroscopic coordinate y G {a,b), converges uniformly. The same scaling leads to a 
nontrivial almost sure limit also for the conditioned motion from [an\ to \ bn\ , i.e. the first 
case of the theorem. The only change is that x is replaced by y on the right hand side of 
the equation ( l43l) . For the paths that are forced to return to their starting point, however, 
no such limit exists because the expectations with respect to pL^nJ^LH depend strongly on 
the transition rate v^anl at the boundary, cf. equation ( l36l) . 

Example. The travel time of a particle through the whole system in a gas with the attractive 
interaction g{k) = I{k > 0} satisfies 

n + T[„] — > / dy 

Jo ^-uy- v{l - y) 

2-U-V (l-M)(l-f), 1-u 

Taking v = and sending m — > yields the value 1/6, which is of course the same as one 
would get by calculating the expected travel time in a gas of noninteracting particles, that 
is g{k) = k. Notice that, somewhat surprisingly, the travel time flH|) is finite even for one 
of the boundary rates tending to 1 (i.e. the radius of convergence of Z). 
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2. Strong disorder: Et;. ^ 



oo 



To prove that the integrals of the functions of the tagged particle position converge for 
in domain of attraction of a-stable variables with a G (0, 1), the functional limit theorems 
of Skorohod 2J] are needed. Of course, there is no hope of getting statements concerning 
almost all environments for nontrivial paths of the tagged particle as in the case of weak 
disorder, but the aim is to determine the distributions of the travel functionals. 

It turns out that theorem 2.7 of reference |24] ideally suits our purposes: The theorem 
concerns the distributional convergence of continuous (in the Skorohod topology) functionals 
of sums of independent and identically distributed random variables under the assumption of 
pointwise convergence. Thus only continuity of the integrals over the macroscopic intervals 
(a, h) must be proven. But this follows from the almost sure continuity of the function 
that maps the paths of Sn^a to the integrand (all this written again in the macroscopic 
coordinates) in the Skorohod topology, and the fact that the Skorohod topology is finer 



than that given by the metric \<t>{x) — '?/'(x)| da; 21|. So the functionals of the tagged 
particle process speeded up by converge as 

^(l + l/a) ^ 

i?L-J ^ f Ql^i+i/. J ^ ^ fiy)Gjix, y) m{Ay) (45) 

in distribution, with Gj{x,y) = {S{x Ay) — S{a)){S{b) — S{y V x))/{S{b) — S{a)) and 
fh^dy) = [Z'{(f){y))/Z{(l){y))]dy. Here S is again the a-stable subordinator and (j){x) = 
uS{x)/ S{1) +f(l — S{x)/S{l)). As in the case of weak disorder, the convergence holds for 
asymmetric paths through intervals as well. The "propagator" must then be substituted 
with Gj{y,y) = {S{y) - S{a)){S{b) - S{y))/{S{b) - S{a)). 

The expectations over the disorder (with respect to P) of the travel time integrals do not 
exist even in finite systems for a < 1/2. To see this, notice that the coefficients Z'(0j)/Z(0j) 
are bounded in the sums to be evaluated, but the expectations of the propagators G j in 
theorem 2 exist for a > 1/2 only: In the simplest case of the set J consisting of just one 
site, 

,,-1 ,,-1 1 /-oo foo -a-l -a-l 

E Y 1 = E / / dxdy 

Vi-i + Vi Vi_i + Vi Ji Ji X ^ + y ^ 



^^<l + x))--T-P^(-^)\x, 
^ " ^ k\{a + k) Vl +a; / 



fc>0 
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with (a;)fc| = a{a + 1) ■ ■ ■ (a + — 1). The sum being uniformly bounded in x, the result 
follows. The lack of convergence of expectations for a < 1/2 is inherited by larger finite 
systems and, eventually, by the continuum limit. Indeed, using the same machinery as for 
calculating the moments and the two-point correlations of the fugacity in section IIIBl we 
get 

EGjix,y) = ar(2 - l/a)(6 - a)^/"-^1>(l)^/" T-^"^^ ' ^'^^^ 

for X < y. Here the gamma function diverges as a — >■ 1/2 from above. 

Example. Let us consider our canonical interacting system with g{k) = I{k > 0} and, 
in particular, the expected time it takes for a particle to travel from the left to the right 
reservoir. Average over the media is also taken. We assume that a G (1/2, 1). Then 

hm n-(^+V«)Ei^o^n+iy =J-Et (Sil) - Siy))J2^' ffmV 

n^oo u-V Jo \S{1)J 

where z = {u — v) /{I — v). As in the proof of theorem 1, we apply Faa di Bruno formula 
in the evaluation the moments of S{y)/S{l). In this case, the resulting summations can be 
carried out explicitly because, instead of a sum over partitions of a fixed number k, we have 
sums over partitions of numbers whose magnitude is determined by the exponential weight 
z^, and the order of these summations can be interchanged. As a result, the expected travel 
time equals 

f\l-y) r«r"-2$(l)e-^"*«^— e-^^'^*««i-^)'^-i)drdy. 
u-v Jo Jo dz 

Performing the derivative and the integrations finally yields the result 
a3$(i)i/-r(2 - l/a)[(l - v){l - uT-' /2 _ ^ _ ^ _ a(u-v)[{l-vr + {l~ur] ^ ^^^^ 



(l-a2)[(l-t;)"-(l-n)"]2 I (i _ _ (i _ ^) 

for it. Taking a ^ 1 leads, up to a model dependent factor $(1), to equation ( |44|) of the 
previous example. Nevertheless, the travel time (1471) diverges for any a G (1/2, 1) as one of 
the boundary rates approaches the critical value 1. 



IV. DISCUSSION 



In the preceding sections, we have shown that the stationary state of a boundary driven 
zero-range process in symmetric random media, with divergent expected inverse jump rates, 
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is in the large system limit described by a normalized a-stable subordinator. It was also 
argued that, in this limit, the physical quantities related to the motion of a suitably speeded 
up tagged particle on an interval J can be calculated using the formula Gj{x,y) = {S{x A 
y) — S{a)){S{b) — S {x y y)) / {S (b) — S{a)) for the propagator or "Green's function", and the 
measure m{dy) = Z'{(j){y))/Z{(f){y)) dy. In mathematics, such function S is known as the 



scale function and the measure m as the speed measure [35|, |36|] of a diffusion. In fact, they 
define a one-dimensional diffusion even when the stochastic differential equation description 
is purely formal, e.g. with discontinuous coefficients as in our case. In practice, the scale 
function and the speed measure tell how the spatial and time axes must be stretched for the 
motion of the particle to be indistinguishable from a Brownian motion (or vice versa, how 
the motion of the particle can be constructed from a Brownian motion). 

Kawazu and Kesten |13i| have proved an invariance principle for a random walk in sym- 
metric random media. The scale function and the speed measure of the motion on large 
scales coincide with S and m given in this article. The speed measure is just the Lebesgue 
measure in this case. In the article of Kawazu and Kesten, the discrete random walk was 
on the infinite set Z, so the only difference to our setting in the noninteracting case is that 
our particles are absorbed at and n -\- 1. For a convergence proof for zero-range processes 
on tori with fixed rates, see |l6|. It would be of interest to find a proof that the speeded up 



motion of a tagged particle in a zero-range process in symmetric random media converges 
to a diffusion defined by the scale function and speed measure given in this article. 



Acknowledgments 

The author would like to thank Stefan Geiss, Pekka Kekalainen, and Juha Merikoski 
for discussions. This work has been supported by the Centre of Excellence program of the 
Academy of Finland. 



20 



[1] F. Spitzer, Interaction of Markov Processes, Adv. Math. 5: 246-290 (1970). 

[2] E. D. Andjel, Invariant measures for the zero-range process, Ann. Prob. 10:525-547 (1982). 
[3] C. Kipnis and C. Landim, Scaling Limits of Interacting Particle Systems (Springer, Berlin, 
1999). 

[4] L. Bertini, A. de Sole, D. Gabrielli, G. Jona-Iasinio, and C. Landim, Macroscopic Fluctuation 

Theory for Non-Equilibrium Stationary States, J. Stat. Phys. 107:635-675 (2002). 
[5] M. R. Evans and T. Hanney, Noncquilibrium statistical mechanics of the zero-range process 

and related models, J. Phys. A: Math. Gen. 38:R195-R240 (2005). 
[6] E. Levine, D. Mukamel, and G. M. Schiitz, Zero-Range Process with open Boundaries, J. Stat. 

Phys. 120:759-778 (2005). 
[7] M. R. Evans, Bose-Einstein condensation in disordered exclusion models and relation to traffic 

flow, Europhys. Lett. 36:13-18 (1996). 
[8] J. Krug and P. A. Ferrari, Phase transitions in driven diffusive systems with random rates, J. 

Phys. A: Math. Gen. 29:L465-L471 (1996). 
[9] S. Havlin and D. Ben-Avraham, Diffusion in disordered media. Adv. Phys. 36:695-798 (1987). 
[10] B. D. Hughes, Random Walks and Random Environments, Volume 2: Random Environments 

(Oxford University Press, 1996). 
[11] O. Zeitouni, Random walks in random environments, J. Phys. A: Math. Gen. 39:R433-R464 

(2006). 

[12] M. B. Isichenko, Percolation, statistical topography, and transport in random media. Rev. 

Mod. Phys. 64:961-1043 (1992). 
[13] K. Kawazu and H. Kesten, On Birth and Death Processes in Symmetric Random Environment, 

J. Stat. Phys. 37:561-576 (1984). 
[14] P. Ferrari, Limit theorems for tagged particles, Markov Process. Related Fields 2:17-40 (1996). 
[15] E. Saada, Processus de zero-range avec particule marquee, Ann. Inst. H. Poincare Probab. 

Statist. 26:5-17 (1990). 

[16] P. Siri, Asymptotic behaviour of a tagged particle in an inhomogeneous zero-range process. 

Stochastic Process. Appl. 77:139-154 (1998). 
[17] S. Sethuraman, On diffusivity of a tagged particle in asymmetric zero-range dynamics, Ann. 



21 



Inst. H. Poincare Probab. Statist. (2006). 
[18] B. Derrida, J. L. Lebowitz, and E. R. Speer, Large Deviation of the Density Profile in the 
Steady State of the Open Symmetric Simple Exclusion Process, J. Stat. Phys. 107:599-634 
(2002). 

[19] K. Sato, Levy Processes and Infinitely Divisible Distributions (Cambridge, 1999). 

[20] I. A. Ibragimov and Yu. V. Linnik, Independent and Stationary Sequences of Random Variables 

(Wolters-Noordhoff, Groningen, 1971). 
[21] P. Billingsley, Convergence of Probability Measures (Wiley, New York, 1968). 
[22] I. I. Gihman and A. V. Skorohod, The Theory of Stochastic Processes I (Springer, Berlin, 

1974). 

[23] A. V. Skorohod, Limit theorems for stochastic processes. Theory Probab. Appl. 1:261-290 
(1956). 

[24] A. V. Skorohod, Limit theorems for stochastic processes with independent increments. Theory 

Probab. Appl. 2:138-171 (1957). 
[25] J. Bertoin, Levy Processes (Cambridge, 1996). 

[26] J. Pitman, Combinatorial Stochastic Processes: Ecole d'Ete de Probabilites de Saint-Flour 
XXXII - 2002 (Springer, Berlin, 2006); J. Pitman, Combinatorial Stochastic Processes, 
Technical Report No. 621, Department of Statistics, University of California, available via 
http://www.stat.berkeley.edu/tech-reports/ 

[27] E. Regazzini, A. Lijoi, and I. Priinster, Distributional results for Means of Normalized Random 
Measures with Independent Increments, Ann. Statist. 31:560-585 (2003). 

[28] J. Gurland, Inversion Formulae for the Distribution of Ratios, Ann. Math. Statist. 19:228-237 
(1948). 

[29] M. Perman, J. Pitman, and M. Yor, Size-biased Sampling of Poisson Point Processes and 
Excursions, Probab. Theory Related Fields 92:21-39 (1992). 

[30] J. Pitman, Some Developments of the Blackwell-MacQueen Urn Scheme, in Statistics, Prob- 
ability and Game Theory; Papers in honor of David Blackwell, edited by T.S. Ferguson, L.S 
Shapley, and J.B. MacQueen (Institute of Mathematical Statistics, Hayward, CA, 1996). 

[31] J. Pitman and M. Yor, The Two-Parameter Poisson-Dirichlet Distribution Derived from a 
Stable Subordinator, Ann. Probab. 25:855-900 (1997). 

[32] J. Pitman, Poisson-Kingman Partitions, Technical Report No. 625, Department of Statistics, 

22 



University of California, available via http://www.stat.berkeley.edu/tech-reports/ 



[33] M. Talagrand, Spin Glasses: A Challenge for Mathematicians (Springer, Berlin, 2003). 

[34] B. Derrida, Random-energy model: An exactly solvable model of disordered systems, Phys. 

Rev. B 24:2613-2626 (1981). 
[35] L. C. G. Rogers and D. Williams, Diffusions, Markov Processes and Martingales, Volume 2: 

ltd Calculus, 2nd Edition (Cambridge, 2000). 
[36] D. Revuz and M. Yor, Continuous Martingales and Brownian Motion, 3rd Edition (Springer, 

Berlin, 1999). 

[37] K. Kook and R. Serfozo, Travel and Sojourn Times In Stochastic Networks, Ann. Appl. Probab. 
3:228-252 (1993). 

[38] F. Baccelli and P. Bremaud, Elements of Queueing Theory, 2nd Edition (Springer, Berlin, 
2003). 

[39] W. Whit, A Review oi L = \W and extensions, Queueing Syst. 9:235-268 (1990). 
[40] R. Serfozo, Introduction to Stochastic Networks (Springer, New York, 1999). 



23 



